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A variation of Rosenstock's trapping model in which TV independent random walkers are all 
, initially placed upon a site of a one-dimensional lattice in the presence of a one-sided random distri- 

bution (with probability c) of absorbing traps is investigated. The probability (survival probability) 
$jv(i) that no random walker is trapped by time t for TV 3> 1 is calculated by using the extended 
Rosenstock approximation. This requires the evaluation of the moments of the number Sjv(i) of 
distinct sites visited in a given direction up to time t by TV independent random walkers. The Rosen- 
i | , stock approximation improves when TV increases, working well in the range Dt In 2 (I — c) <C In TV, D 

• being the diffusion constant. The moments of the time (lifetime) before any trapping event occurs 

are calculated asymptotically, too. The agreement with numerical results is excellent. 
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I. INTRODUCTION 



Survival of Brownian particles in a medium populated with randomly distributed static traps is a fundamental 
problem (the "trapping" problem) of random walk theory that has been an active area of research for decades with 
many applications in physics and chemistry The origin of this problem can be traced back to Smoluchowski's 

theory of coagulation of colloidal particles flll-H]. It has now become a basic model of widespread interest in areas 
such as trapping of mobile defects in crystals with point sinks [^,3|, the kinetics of luminescent organic materials ||, 
the kinetics of photosynthetic light energy to oxygen conversion [7], anchoring of polymers by chemically active sites 
H , atomic diffusion in glasslike materials || and many more |l(J . 
This paper is devoted to a variation of the so-called Rosenstock trapping problem on a one-dimensional substrate, 
ly-j . Usually, the one-dimensional Rosenstock trapping problem (which we will call the "two-sided" Rosenstock trapping 
' problem for reasons that will be apparent later on) is stated as follows |j]-|3| . A one-dimensional lattice is filled with a 
random distribution of static traps; then, one (TV = 1) random walker is placed initially (t = 0) at a given site of the 
lattice; it starts to diffuse and, eventually, is caught by a trap. In this paper we study a different but closely related 
trapping problem (which we will call the "one-sided" Rosenstock trapping problem) in which (i) only a half-line of 
a one-dimensional lattice is filled with a random distribution of static traps with concentration c (this process could 
mimic the excitation or production of defects in one side of a fiber by irradiation, the other side being shielded) and 
(ii) TV independent random walkers are placed initially (t = 0) at the contact point (x — 0) between the two half-lines 
that is taken as origin. These random walkers start to diffuse, and eventually one of them is trapped at the nearest 
site occupied by a trap (or deactivates it). The statistical quantities of main interest are the survival probability, 
$Ar(i), defined as the probability that no random walker has been trapped by time t, and the lifetime, T/v, defined as 
the average time at which the first random walker of the set of TV arrives at a trap site. To our knowledge, this is the 
first multiparticle (TV ^ 1) Rosenstock trapping problem ever studied. A good reason for this is that the multiparticle 
versions of the trapping problem are much more difficult to solve than the trapping problems with a single particle. 
This fact will be evident is this paper: the present trapping problem is elementary for TV = 1 and, in this case, we 



will report the main results for the sake of completeness only (see section III). On the other hand, the multiparticle 
version is much more involved (see section |^ and 0). As an exact evaluation of 3?jv(i) for TV > 1 is elusive, we 
have resorted to asymptotic analysis techniques. In particular, for TV ^ 1, we have used the extended Rosenstock 
approximation (or truncated cumulant approximation) [|l|J^,|lQ| p~2| - This requires to find the first moments (S^}(t)), 
m = 1,2,.. ., of S^it), the number of sites situated to the 'right' of the origin x — that were visited up to time t 
by TV random walkers which started at the site x — at time t — 0. Note that Sjv(£) is just the maximum distance 
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reached by any of the N random walkers in the +x direction from the origin by time t, i.e., S/v(i) is the one-sided 
span of the iV-particle random walk. The problem of evaluating the first moment (SV(i)) has already been addressed 
in |l3|-[l5| explicitly and in 16 1S|] implicitly. The quantity studied in these last references was the number of distinct 
sites explored by N random walkers on the one-dimensional lattice in either direction, SW(i), but (iSjy( t)) = (5 l jy(t))/2. 
However, little is known about higher moments of Sjv(i), (S™(t)), except for some rough estimates . The 

idea of evaluating the survival probability for the multiparticle trapping problem by using the moments of Sn(£) into 
the Rosentock approximation was suggested by Larralde et al. in jl7|, although, to the best of our knowledge, it has 
not been implemented perhaps for the lack of precise expressions for the moments (S™(t)). Therefore, in order to 
use the Rosenstock approximation in our one-sided multiparticle trapping problem, we must find rigorous asymptotic 
results for the moments of SN(t). This is another objective (important in itself) of this present work that, besides, 
we hope can illuminate how to deal with the evaluation of the moments of Sjsr(t) for other substrates. 

The paper is organized as follows. In section || we recall the connection between the survival probability $Ar(i) 
and the so-called fixed trap survival probability (i.e., the survival probability when the trap is placed at a given 
distance) and give the basics of the extended Rosenstock approximation. The one-sided trapping problem for N = 1 
is addressed in section H. In section IV we calculate the moments of 5;v(t) in the form of an asymptotic series in 
which the corrective terms decay as powers of 1/lniV. The evaluation of (Sp}(t)) is a necessary prerequisite for the 
implementation, in section V| of the Rosenstock approximation for <&jv(^)- In this section we compare $jv(i) as given 
by the extended Rosenstock approximation with numerical results. The moments of the lifetime Tn and its variance 
are evaluated in section W\, The paper ends with some conclusions and remarks. 



II. ONE-SIDED TRAPPING MODEL AND ROSENSTOCK APPROXIMATION 



The one-sided Rosenstock trapping model is defined as follows: (i) Quenched traps are randomly distributed on 
the right-hand side of a one-dimensional lattice (x > 0) with concentration c (1 — c = p); (ii) the random walkers 
are placed initially upon site x — which divides the randomly filled trapping half-line and the empty one; and (iii) 
the traps are irreversible, that is, a walker encountering a trap is killed there. Then, the survival probability is given 
by $jv(*) = Y^r^iP 1 PnIA 7 ") j where P N (t\r) is the probability that the span of the N random walkers in the positive 
direction (the largest distance reached by any of the N random walkers for x > 0) is equal to r after t time steps [0 . 
Let Tjv(t|r) be the probability that the site x = r has not been visited by any of the N random walkers by time t 
(the so-called fixed-trap survival probability). Then, in the continuous limit, 



' r dT N {t\r) j 

p : dr 

dr 



(1) 



where we have used the relationship Ppf(t\r) — dT ^{t\r) / dr between the one-sided span distribution P^{t\r) and 
the fixed-trap survival probability r^v(t|r). For N independent random walkers one has rV(i|r) = [T(t|r)] where 
r(t|r) = ri(i|r) is the probability that distance r has not been reached by a single random walker by time t. For the 
one-sided diffusion process it is well-known that [|],|| : 



r(t|r) 



(2) 



The extended Rosenstock approximation (or truncated cumulant expansion) is now a standard approach jl] [| to 
the Rosenstock trapping problem, which was first proposed by Zumofen and Blumen |Tl| ] and that we recall here for 
the sake of reference. From the definition of S^ft) as the number of distinct sites on the positive half- line visited up 
to time t by N independent random walkers that started at x = at time t = (note that this means Sjv(0) = 0), the 
survival probability of the N random walkers is $jv(i) = (p SN ( 1 ^ = (exp[SV(i) hip]) . The average in this equation 
is performed over all realizations of the random walkers' exploration of the lattice up to time step t. The well-known 
cumulant expansion technique [00] allows an alternative form of $Ar(i) as an infinite series expansion 



$jv(t) = exp 



K n (lnp) r 

71=1 



n 



(3) 



where K n , n = 1, 2, . . . denote the cumulants of Sw(t): Ki — (5jv(i))> K 2 = (Sjy-(i)) — (SN(t)) 2 = o%{t), ... If only the 
first term of the sum in Eq. ([$]) is kept we arrive at the zeroth-order Rosenstock approximation 



$W(£) = e <S«(*)>m P . 



(4) 
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The error made by using this approximation can be estimated by taking into consideration the next exponential 

term in Eq. (§): ® N (t) = (t) 1 + O (o% (Inp) 2 ) . Thus, the condition a 2 N < l/(lnp) 2 must be fulfilled for the 

zeroth-order Rosenstock approximation to be reasonable. The first-order (extended) Rosenstock approximation is 
obtained by retaining two terms in the infinite sum of Eq. (||) : 



3$(t)=exp 



(S N (t))\np+-a 2 N (lnp) 2 



(5) 



The relative error of this expression is of order O [^(lnp) 3 ] . 



III. ONE-SIDED TRAPPING PROBLEM WITH A SINGLE RANDOM WALKER 



The one-sided trapping problem is quite simple for N = 1 . Anyway, we report here the main results for the sake 
of completeness. From Eqs. (TP) and (0) the survival probability of a single random walker is 



drp r —erf 
dr 



tL= ) = e ^ erf c ( X - 



(6) 



where x = ^/ADt ln(l/p). For very long times, x — > oo, the asymptotic expansion of the complementary error function 
[Il9l allows us to write 



(2m- 1)!! 



Jim 



where (2m — 1)!! = (2m — 1) ■ • -5 • 3 • 1. Thus, an asymptotic time regime is reached for t S> l/(lnp) 2 where 
the survival probability exhibits a power law decay <&i(t) « l/[y/nD ln(l/p)]i -1 / 2 . This is an algebraic fluctuation 
slowdown corresponding to the Donsker-Varadhan limit. 

In order to apply the extended Rosenstock approximation to the single random walker case for small x, we must 
evaluate the moments of the one-sided span Si (t) : 



(8) 



or, after integrating by parts, 

(Sr(i))=m 



erfc 

o W^DtJ 



r m - L dr 



m + 1 



(ADt) 



m/2 



(9) 



with r(m) being the gamma or factorial function and where Eq. (g) has been used. It is now easy to verify that 
n n = a n (4Dt) n / 2 , n= 1,2,... with a x = l/y/n, a 2 = 1/2- 1/tt, a 3 = (4 - n) / (2w^/tv) , a 4 = (2tt - 6)/tt 2 , .... Direct 
substitution of these cumulants into the general expression for the Rosenstock approximation (|^) yields 



= exp ]T ^ m) n/2 ^ n P>= exp £ 



(-1)" 



exp 



1 



(10) 



7T-2 



7T - 4 



,n=l 

7T - 3 



.x 4 + G(x 5 ) 
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Notice that the Rosenstock approximation given by Eq. (jl(]) coincides with the exact result in Eq. 



IV. ONE-SIDED SPAN OF A SET OF RANDOM WALKERS 



The objective of this section is twofold: first we want to obtain rigorous asymptotic expansions for the one-sided 
span moments 



3 



dT N (t\r) 
dr 



(11) 



for m = 1, 2, . . . and N ^> 1 independent random walkers, and, second, we want to check the reliability of the obtained 
asymptotic expressions comparing them with numerical results. 
Integrating Eq. (fijl) by parts, we find 



/>oo 

(S£(t)) = m^Di) 1 "/ 2 / r _1 [1 - erf"®] d£. 
Jo 



(12) 



where £ = r/v4~Dt. In order to evaluate this integral for large values of N it suffices to know T(t|r) = erf(£) for large 
£, namely, erf(£) » 1 — 7r _1 / 2 £ _:1, e _ £ (l — £~ 2 /2 + •■•). The asymptotic evaluation for large N of the integral of Eq. 
( |l2| ) is not a easy task. Fortunately, if one compares this integral with the one carried out in Ref. fL8|| , one realizes 
that both integrals are formally equivalent. In this way one finds 



(S%) « [4^1n(7V)] m/2 [l-A(m)] 



(13) 



where 



A(m) = Y J ln~" N ^ s<f } (m) ln J In N. 



(14) 



n=l 



Up to second order (n = 2) the coefficients are 

(l)/ \ 
s >(m) = - 

s-[ (m) 
4 2) (m) 



muo 
'~2~' 



77! 

7 

77! 

8~ 



(2-m) ( — +cc 2 



2w + 2 



»f ) (m) = ^[(m-2)w-l] 



4 2) M = §(2-m), 



(15) 
(16) 

(17) 
(18) 
(19) 



where u; = 7 — iln7T = 0.0048507 • • • and 7 = 0.577215 • ■ • is Euler's constant. In particular, the average value of 
Sjv(t) i s given by 



(S N ) « [4£>iln(A0] 



1/2 



lnlniV -2w i In 2 In AT - (1 + w) In In N + 5^(1) 



s( 2 )/ 



41niV 



8 In 2 AT 



(20) 



with §q (1) = 7r 2 /6 + o> 2 + 2cj + 2 = 3.654659 • • •. It is remarkable how close the result obtained in ||15[, namely, 



{S N )^[4Dt\n(N)} 1/2 



1 - 



lnlniV + ln(4?r) 
4 In AT 



(21) 



is to the rigorous result given by Eq. (^0|) . Notice that the result of Eq. ([H]) starts to differ from that of Eq. (p0| ) in 
the first corrective term: in this equation ln(47r) = 2.531024 • • • plays the role of — 2uj = —0.009701 • • • in Eq. (EOf). 
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FIG. 1. The dependence on N of the first seven moments of Sjv(t). We have plotted S = {S^} /[4Dt In N] m/ 2 versus 1/ln/Y. 
The symbols correspond to the numerical estimate for TV = 2" with n = 3, 5, . . . 30, when m = 1 (open circles), m = 2 (filled 
circles), m = 3 (open squares), m = 4 (filled squares), m = 5 (open triangles), m — 6 (filled triangles) and m = 7 (crosses). The 
solid lines correspond to the second-order asymptotic approximation 1 — A(m) as given by Eq. (114). Notice the importance of 
the corrective terms as the main term (zeroth-order approximation) would be an horizontal line passing by S = 1 
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FIG. 2. The ratio R(m) = [(Sn) /o"jv(m)] m / 2 versus In N for (from top to bottom) m = 2,3, ... ,7. The symbols correspond 
to the numerical evaluation for N = 2 n with n = 0, 1, . . . 30. The lines correspond to the asymptotic result of Eq. ([^J). The 
slope of the linear fit to the last five points (N from 2 26 to 2 30 ) is 1.53 (1.56) for m = 2, 0.88 (0.90) for m = 3, 0.62 (0.64) for 



m = 4, 0.48 (0.49) for m = 
asymptotic value {48/ [m(m 



5, 0.39 (0.40) for m = 6, and 0.33 (0.34) for m = 7, in good agreement with the corresponding 
- l)^ 2 ]} 1 ' 2 given by Eq. (|23j) which we have set in parentheses. 



In Fig. |]we compare the values of {S 1 ^) /[4J5t In TV]™ 1 / 2 for m = 1, 2, ... 7 obtained by integrating Eq. ( fl2"l) numerically 



with those obtained by means of the second-order asymptotic approximation of Eq. (Uq). The importance of the 
corrective terms given by A is evident as well as the good performance (even for not-too-large values of N) of the 
second-order asymptotic expression. This is especially notable for low-order moments. Notice that, at least for 
n = 1, 2, the coefficient s^\m) is a polynomial of degree n in m, and that, if this property holds for all n, it would 



explain why the truncated asymptotic expansion of (S^) 
From Eq. (|l|) one finds 



worsens for increasing values of m. 



a N (m) = (S%) - (S N ) 



48 



(m - l)7r 2 (4Di) m / 2 (ln/V) m / 2 - 2 



l + O 



In 6 In AT 
In TV 



(22) 



(TAr(m) m(m — 1)7T 2 



(S. 



48 In 2 N 



l + O 



ln J In N 
In TV 



(23) 
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Figure ^ is a plot of the numerical and analytical results for R{m) = [(Sn)" 1 /crN('m)] 1 ^ 2 versus In AT. It is clear 
that the numerical results closely follow the theoretical prediction given by Eq. ( p3| ) although some differences are 
noticeable. We attribute the difference to the existence of corrective terms of order In 3 In N/ In N . 
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FIG. 3. The ratio R(2) = (v24/7r) ln(iV) (7ff/{Sjf) versus 1/lniV. The symbols correspond to the numerical evaluation for 
N = 2 4 , 2 5 , . . . , 2 30 . The line is a curve of the form a(l + ln~ 1 N ^] ,_ Q bj In-' In N) fitted to the last fifteen simulation points, i.e., 



the points corresponding to N — 2 
and 6 3 = 0.005. 



The fitting parameters are a = 0.998, b = 0.009, 61 = -0.019, b 2 = 0.094 



Let us look at this question more closely. In Fig. we have plotted the ratio 



R(2) 



24 H N)- aN 



(S. 



l + O 



N 



ln d lniV 
In TV 



(24) 



evaluated numerically versus 1/lniV, where cr^ = <tjv(2) is the variance. We note that, at first sight, there are two 
features that cause unease in this figure: first, the numerical results are relatively far from unity (the main asymptotic 
term) even for very large values of N, and second, it is difficult to state that unity is the final value for N — > 00 by 
only looking at the points (the numerical results). We can shed some light on these two aspects by considering the 
form and value of the corrective terms of Eq. ( pl| ) . The first aspect can be understood by taking into account that the 
functions ln n (lniV)/ In N have non- negligible values and are only (slowly) decreasing functions for very large values 
of N (e.g., when n = 3, this function is only decreasing for N > 5 x 10 8 ). In order to resolve the second difficulty, 

we use our knowledge of the form of the first corrective term of Eq. ( p4| ) , namely, Y^=o a< f^ ^ hi(iV) / In N, to fit the 
numerical results of the quantity v / 24(ln N)a^ / (tt(Sn)) to the expression a(l + In -1 NJ2j=o bj m ^ m N). Neglecting, 
for example, the values corresponding to N = 2°, . . . , 2 10 (recall that our expressions are asymptotic expressions valid 
for large N), the fitted function leads to a ~ 0.988, a result that is in excellent agreement with our predicted value of 
unity [i.e., the main term of Eq. (p4|)]. 

Finally, it is interesting to note that Eq. (p2|) tells us that, up to first order in 1/lniV, we can get the TO-th moment 



(Spf) from only the knowledge of the first moment (SW) 
through the relationship 



To be more precise, let us define the parameters sj™' (m) 



(S. 



[4Dtln(N)] m/2 [l-A(m)] 



(25) 



where 



A(m) = ^2 ln 



~ n N Vsl n) (m) ln J ln7V. 



3=0 



(26) 



Hence, using Eqs. (p^|)-(p^|), we find (i) that s^ l \m) 
and (ii) that the first terms of (S^) and (5 



to j = and 



.-(2) 



(to) 



N/ 
.(2) 



(m) = to(to 



Sj(m) for j — 0, 1 when n = 1 and for j — 1,2 when n = 2, 
that are different are second-order terms, namely, those corresponding 
- 1)tt 2 /48. 



It is also worthwhile noting that, working up to second-order asymptotic corrective terms (n = 2), the j-th cumulant 
with j > 2 is zero, or, in other words, that the distribution of Sn is, up to this second asymptotic order, Gaussian. 
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V. ONE-SIDED TRAPPING PROBLEM WITH A SET OF RANDOM WALKERS 



From Eq. (Jl]) and after integrating by parts, we can write the survival probability of a set of N independent random 
walkers in the one-dimensional one-sided problem as 



*jv(t) 



[erf (£)f % 



where x = \J ADt ln(l/p) and £ = r/VADt. 

The asymptotic behaviour of <$>isr(t) for an arbitrary number of particles and x 
Watson's lemma [^0| and the expansion of the error function erf(£) for small £: 



(27) 



oo is a direct consequence of 



N 



Nl 



(N + 2)\N 
3x N + 2 



(28) 



For TV 
$j\r(t) 



= 1 we recover the result in Eq. (Q). An slow power-law time decay of $N(t) is observed in this limit, 



t- N ' 2 . 

It is interesting to note that the present one-sided multiparticle Rosenstock trapping problem is related to one of 
the predator-prey problems discussed by Krapivsky and Redner |l4| , |l5f| in which a static prey or "lamb" is captured 
by one of a set of N diffusing predators or "pride of lions" . These authors considered the case of TV predators and 
a prey at given relative positions so that the case we study here differs in the sense that the traps (or preys) are 
randomly distributed. In their analysis, Krapivsky and Redner found that the survival probability of the "lamb" is 
given by the power law t~ N / 2 . So, this behaviour agrees with that of our stochastic "lamb" problem in the long-time 
regime, x 3> 1, i.e., t ^S> lnTV/[ln(l — c)] 2 . The reason for this behaviour is that, for very long times, how the traps 
are configured is not essential in the trapping kinetics, and the the slow $jv(i) ~ t~ N / 2 power law decay of the fixed 
tramp (or "lamb") case settles down. 

Next, we will deal with the asymptotic behaviour of $Ar(i) for small x and large TV. In order to implement 
the Rosenstock approximation, we will use the expressions for the average one-sided span (Sjv(i)) and its variance 
a N (2) = a 2 N 



given by Eqs. (^0|) and (|22|). Then, the zeroth-order Rosenstock approximation given by Eq. ([|) is 



I - 



In In AT - 2lo 



4 In AT 



(4DilnTV) 1/2 ln(p) 

\ In 2 hxN - (1 + w) In In N + tt 2 /6 + u? + 2lo + 
8 In 2 



2 /hi 5 In AT 
ln J N 



(29) 



In section || it was mentioned that the relative error for is of order O (a% ln 2 p), so that the condition a 2 N <C 
l/(lnp) 2 must be fulfilled for the zeroth-order Rosenstock approximation to be reasonable. Using the expression for 
<j 2 N given by Eq. (p2h, one sees that the relative error goes as 7r 2 Dt/[6(lnp) 2 In N], so that 4>^(i) should give good 



results when t <C lnA^/(lnp) 2 or x <C Vhi N. Consequently, there exists a time regime, which becomes larger as 
the number of random walkers increases and the concentration of traps decreases, where the approximation for the 
survival probability in Eq. ( |29| ) must work reasonably well. We will denote by ^^ m \t), m = 0, 1, . . . , the zeroth-order 
Rosenstock approximation for the survival probability that results from retaining m corrective terms in the expression 
for the average one-sided span and ignoring the rest of the terms. For example, 



$^ 1) (i) = exp 



- a; (In AT) 



1/2 



f 



In In AT - 2lu 
4 In TV 



(30) 



Keeping all the explicit terms in Eq. (p9|), we get <j>^(f). 



Proceeding as above, the approximation given by Eq. (^) that includes the variance term, i.e., the first-order 
Rosenstock approximation, is 



In = (4DtlnA^) 1/2 ln(p) 



f 



lnlnAT - 2lu 



41n Af 

In 2 In N - (1 + uj)\nlnN + tt 2 /6 + lu 2 + 2lu + 2 



I 1-2 

4 



(4DtlnAT) ln 2 (p) 



8 In 2 N 



48 In N 



O 



hfhiN 
In 3 N 



(31) 



O 



\n A In TV 
In 3 A^ 
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We have seen in Sec. || that one can estimate the error of the hrst-order Rosenstock approximation by looking at the 
value of ^(lnp) 3 . The kurtosis, K3, was not calculated in section IV explicitly, but we know (see the last paragraph of 

that section) that it is zero up to at least second-order corrective terms, i.e., at least, K3 = O ft 3 / 2 In 3 hi7Vln- 3/2 TvY 

This means that the relative error of the first-order generalized Rosenstock approximation is a quantity of order 
0[(x/\/\n TV) 3 In 3 In TV] which will be much smaller than that corresponding to Eq. (||), i.e., 0(x/y/hi TV), as long as 
the condition x <C Vhi TV is fulfilled. 




15.0 



FIG. 4. Survival probability of TV = 10, TV = 10 3 and N = 10 6 random walkers in the one-dimensional one-sided Rosenstock 
trapping model versus x — V '4-Diln(l j 'p) . The circles correspond to a numerical integration of the exact result, the broken 



lines correspond to the Rosenstock approximations of orders > 



<1> 



(01) 



(02) 



enumerated from below, and the solid line 



corresponds to $^' ) . The lines corresponding to Q^ 1 ' and for TV = 10 5 are indistinguishable 



(01) 



--,(02) 



In Fig. ^ we compare the different order Rosenstock approximations with the numerical evaluation of Eq. ( p7| ) for 
TV = 10, 10 3 and 10 6 . We observe that adding more corrective terms increases the agreement of the Rosenstock 
approximation with the numerical result for small values of x. Quite noticeably, the approximations also get better 
as TV increases. 



VI. LIFETIME IN THE ONE-SIDED TRAPPING PROBLEM 



Let the lifetime of the set of TV independent random walkers T/v be defined as the time at which some random 
walker of this set is first trapped or, conversely, the time at which the lamb is killed by the pride of TV lions if the 
expression coined by Krapivsky and Redner is used p^ , |l5| |). The m-th moment of the lifetime distribution is given 
by 



(T 



N , 



or, taking into account Eq. (p7|), 



dt 



(32) 



ro=m(l/» / drp r t™(r) 
Jo 



(33) 



where t^(r) is the m-th moment of the time to first reach the distance r by the first random walker of a set TV 
independent diffusing random walkers: 



C JVV 



•)=- 



00 d 
t m — 
dt 



N 



dt 



(34) 



Of course, t^(r) can be understood as the moments of the lifetime in the trapping problem with a fixed trap at 
a distance r from the starting site of TV independent random walkers. This problem has been widely studied. For 
example, Lindenderg et al. pl| studied the first passage time for small TV, finding that the m-th moment of the first- 
passage-time distribution for the first of the walkers to reach r is finite if there are at least 2m + 1 random walkers 
starting from the origin. Hence, by using Eq. ([33]), we can extrapolate these conclusions to the one-dimensional 
one-sided Rosenstock's trapping problem. In particular, this means that (Ti) and (T2) are infinite but (Tjy) = 



Cjv/(2-Dln 2 p) is finite for every TV > 3, the coefficient Cn being given by J °° [erf(£ 



- 1 / 2 )} N dt. 
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The trapping problem with fixed trap and large N has also been studied ]22| . |23| . In particular, for the one- 
dimensional lattice 



4D 



1 



1 



where Aq = Then, from Eqs. 



\n m X N t InAoiV 
and ((351) we have 



1 



lnlnAoiV - 7 + 



(T$) = t 



N 



( [(2m)!] 1 



/2r. 



\ In 1/p 



(35) 



(36) 



Thus, the m-th moment of the lifetime {T^) coincides with the m-th moment of the first passage time to a fixed trap 
at a distance [(2m)!] 1//2m times the average distance (I) = — 1/lnp of a trap to the origin in the one-sided Rosenstock 



model with trap density c — 1 — p. According to Eqs. (36) and (p5|), and writing the expansion in terms of In TV 
instead of In XqN, we finally have 



(2m)! 



[4D (hip) 2 In TV] 1 



1 



In TV 



- In In AT - uj 
2 



2 In N 



1 



(TO + l)7r 2 



+(m + 1)uj 2 + ui — ( muj + oj + 



In In N 



(m + 1) ln 2 ln/V 



6 



O 



In 3 In/Y 
In 3 TV 



(37) 



The asymptotic expansion for the variance of the distribution of first passage times, a 2 N = (Tjy) — (T/y) 2 , is calculated 
from Eq. @: 



20 



AD (hip) 2 InA^ 



In TV 



- hi In TV - oj] H \ — 

2 J 21n 2 /V 



Stt 2 
15 



+ ( 1 + 2.„- ! - ( - + 2 w J In In N + i In 2 In TV 



(9 



In 3 IniV 
In 3 TV 



(38) 



The large size of the variance is remarkable: notice that the coefficient of variation (Tn)/ctn ~ l/Va (for AT ;g> 1) is 
less than unity. 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 
InN 



FIG. 5. Inverse of the lifetime of a set of TV independent random walkers in the one-dimensional one-sided Rosenstock 
trapping model with a density of traps c = 1/2 versus In TV. Numerical integration results are plotted as circles and the lines 
are theoretical asymptotic predictions of zeroth (dotted), first (broken), and second orders (solid). The value D — 1/2 has been 
used. 



9 



3.0 




2 4 6 8 10 12 14 
IniV 



FIG. 6. The same as Fig. |5| but for the inverse of the square-root of the variance, l/crjv. 

In Figs. || and ^ we have plotted the inverse of the average first passage time, l/(T/v), and the inverse of the 
square-root of the variance, 1/erjy, respectively, versus IniV. The second-order approximations for the first passage 
time in Eq. ( |37| ) and the variance in Eq. ( p8| ) are in excellent agreement with the results of the numerical integration 
ofEq. (§f). 



VII. CONCLUSIONS AND REMARKS 



In this paper we have studied the one-dimensional one-sided multiparticle Rosenstock trapping model in which N 
independent random walkers start their random exploration at the same site, x — 0, of a one-dimensional lattice whose 
right side, x > 0, has been randomly filled with static traps with a density c= 1 — p. This problem can be seen either 
as a multiparticle version of the usual "two-sided" trapping model, or as "random" version of the "predator-prey" 
problems such as that discussed in Refs. fl4|,|l5| regarding the survival probability of the prey, or in Ref. [^lj regarding 
the prey's lifetime. Our main interest has been the calculation of the survival probability <5>N(t) and the moments of 
the lifetime (Tff). 

For evaluating $jv(i) for large N we resorted to the (extended) Rosenstock approximation. This required, as an 
intermediate step but interesting in itself, the asymptotic evaluation of the moments of the one-sided span Sjv(i) 



in terms of a series of the form (S%) » {ADt ln7V) m/2 1 - £,?=i( ln N)- n YT j= a s " lnJ ln N ■ The importance of 
evaluating the corrective terms even for very large values of N is evident as these terms decay mildly as powers of 
l/lniV. It is worth mentioning that we have found that (S™(t)) = (SN(t)) m up to first- order corrective terms (n = 1). 
The agreement of the Rosenstock approximation with numerical results for $jv(0 improves when N increases and 
x = — v4Dtln(l — c) decreases, the agreement being very good when x <C Vln N. However, the stretched exponential 
behaviour of <frjv(i) characteristic of the Rosenstock approximation breaks down for x \/ln N and an algebraic 
long-time tail behaviour settles in: $Ar(i) ~ t~ N/>2 . This is the fluctuation slowdown reported in the trapping problem 
literature p4| and known as the Donsker-Varadhan limit, but here more dramatic, being algebraic instead of stretched 
exponential, because of the infinite half-line free of traps that exists in the one-dimensional one-sided trapping model. 

Once the survival probability was determined, we dealt with the problem of the first passage time TV of the first 
random walker of a set of N to the nearest site occupied by a trap. We found that the m-th moment (T$ ) coincides 
with the m-th moment of the first passage time to a single fixed trap placed at the distance r m = [(2m)!] 1 ^ 2m '(I) from 
the origin, (I) = 1/ ln(l/p) being the average distance to the origin of the nearest trap of the random distribution. 

Finally, it must be remarked that the extended Rosenstock approximation could also be used for the multiparticle 
trapping problem on the "two-sided" one-dimensional lattice and on other Euclidean and fractal substrates because 
asymptotic series for the average number of distinct sites visited (Sjv(i)) nave previously been derived for these media 
pSf . However, the calculation of higher order moments of poses a problem of completely different order of 

magnitude that still remains unsolved. On the other hand, we have found that, up to first-order corrective terms, 
(jSjy-(i)) = (S]y{t)) 2 - At this point, it is tempting to conjecture that this holds for 5jv(t) too. Were this true, the 
variance of Sj\r(t) would be a quantity of second order and, therefore, the zeroth-order Rosenstock approximation with 
(Sn) given up to first-order corrective term should lead to a good approximation for the survival probability $jv(^) 
for Euclidean and fractal media too. Work is in progress to check these conjectures and thereby extend the results of 
this present work to more general trapping problems. 
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